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Abstract 



We investigate the coherent dynamics of carriers in semiconductor su- 
perlattices driven by ac-dc electric fields. We solve numerically the time- 
dependent effective-mass equation for the envelope function. We find that 
carriers undergo Rabi oscillations when the driving frequency is close to the 
separation between minibands. 
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Recent advances in laser technology make possible to drive semiconductor nanostruc- 
tures with intense coherent ac-dc fields. This opens new research fields in time-dependent 



ultra-high speed devices. Artificial two-level semiconductor nanostructures acting as switch- 
ing devices based on Rabi oscillations (RO's) have been already suggested [[J. Moreover, 
resonant phonon-assisted tunneling through a double quantum dot could be used as an 
efficient electron pump from spatial RO's [{§]. 

Zhao et al. have investigated analytically a tight-binding model of a two-band system in 
a time-dependent ac-dc field in the weak coupling limit ||. They identified RO's between 
Bloch bands under resonant conditions, which reveal the existence of quasienergy bands and 
fractional Wannier-Stark ladders. The advances achieved in molecular beam epitaxy, which 
allow to fabricate semiconductor superlattices (SL's) tailored with the desired conduction- 
band profiles, make these systems ideal candidates to propose experiments on coherent 
carrier dynamics. However, the tight-binding approximation presents some limitations to 
describe actual semiconductor superlattices (SL's) when the coupling between neighboring 
quantum wells is not weak. Thus, in order to experimentally access the validity of theoretical 
predictions, one should use a more realistic model. In this Letter we present an effective- 
mass model beyond the tight-binding approximation, containing all ingredients of actual 
SL's, namely finite interband coupling and multiband scattering. 

We consider electron states close to the conduction-band edge and use the effective-mass 
approximation. The electron wave packet satisfies the following equation 



where x is the coordinate in the growth direction, F and u are the strength and the frequency 
of the ac field. The SL potential at flat band is Vsl(x) = AE C if x lies inside the barriers 
and zero otherwise, AE C being the conduction-band offset. We have considered a constant 
effective-mass m* for simplicity. 



transport in mesoscopic systems []T]|| and puts forward the basis for a new generation of 



dt 
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The band structure at flat band is computed by using a finite-element method ||. The 
eigenstate j of the band i with eigenenergy is denoted as ip^{x). A good choice for the 
initial wave packet is provided by using a linear combination of the eigenstates belonging 
to the first miniband. For the sake of clarity we have selected as the initial wave packet 
ty(x,Q) = ip*f\x), although we have checked that this assumption can be dropped without 
changing our conclusions. The subsequent time evolution of the wave packet *f?(x, t) is 
calculated numerically by means of an implicit integration schema 0. In addition to *&(x, t) 
we also compute the probability of finding an electron, initially in the state *&(x, 0) = ip\ {x), 
in the state ip^\x) 

P$(t)= r dx**(x,t)^\x). (2) 

J — oo 

We present here the results for a SL with 10 periods of 100 A GaAs and 50 A Gao.7Alo.3As 
with band offset AE C = 250 meV and m* = 0.067m, m being the free electron mass. We 
consider electric field F = 25kV/cm as a typical value, although similar results are observed 
for other strengths. Figure |I|(a) displays Pqi (t) for F = 25kV/cm at the resonant frequency 
ct; ac = (E[ 5 ^ — Eq 5 ^)/?i = 150 THz. Thus, we are monitoring the transitions between the 
central state (j = 5) in the first miniband to the central state in the second miniband as a 
function of time. We observe the occurrence of very well defined RO's with an amplitude 
close to 0.3. Summing up the probabilities of the rest of states in the second miniband, 
the probability of finding the electron in this band is very close to unity (~ 0.99). The 
frequency of the RO's, obtained performing the fast Fourier transform (FFT) of Pqi (t), is 
^Rabi = 19.18 THz. The probability Poi\t) is dramatically reduced when the ac driving field 
is out of resonance, as shown in Fig. |l|(b) for a; ac = 100 THz. The FFT of those date reveals 
no specific features besides the peak at the driving frequency oj ac . 

In a pure two-level system, a straightforward perturbation calculation yields WR a bi = 
I Foi\/h at resonance, where Fqi is the matrix element of the perturbation between the 
ground state and the first excited state. Thus, i^R a bi is linear in the electric field in a pure 
two-level system. Although the SL is not a pure two-level system, we realize that this linear 
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dependence still holds, as we can see in Fig. 0. 



As an estimation of the leakage current that one could observe in electron pumping 
devices based in RO's in SL's, we have studied the integrated density in the right part of 
the SL, defined as 



where x r is the coordinate of the right edge of the SL. Figure |3| shows the results at the 
resonant frequency uj. ac = 150 THz as well as out of resonance when cu ac = 200 THz. Under 
resonant conditions the wave function is emitted by bursts from the SL region every time a 
RO has been completed. On the contrary, the tunneling across the whole SL is negligible 
when the ac driving frequency is not very close to the resonant one. 

In summary, we have studied the coherent carrier dynamic in semiconductor SL's driven 
by an intense ac-dc field. We found that the electron can perform RO's under resonant 
conditions as in pure two-level systems. The Rabi frequency of the oscillations depends 
linearly on the strength of the electric field. We have shown that electrons are emitted by 
bursts under resonant conditions, whereas the tunneling probability in vanishingly small 
out of resonance. Therefore, we suggest that semiconductor SL's driven by an intense ac- 
dc field may be used as an efficient electron pumping device in THz science. Finally, we 
should mention that we have not considered the role of imperfections or any other scattering 
mechanism that could result in a reduction of the carrier coherence. Our preliminary results 
show that scattering by interface roughness does not strongly modify the above picture. 
Further work along these lines is currently in progress. 
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FIGURES 
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FIG. 1. The probability of finding an electron, initially situated in the state ipQ , in the state 
V'j 5 ^ as a function of time for F = 25 kV/cm, when the ac field (a) is tuned to the resonant frequency 
u; ac = 150 THz and (b) is out of resonance w ac = 100 THz. Note the different vertical scales. 

FIG. 2. Rabi frequency as a function of the dc field when the driving frequency is tuned to the 
resonant frequency w ac = 150 THz 

FIG. 3. Pr{t) as function of time at at the resonant frequency w ac = 150 THz (solid line) and 
out of resonance uj ac = 200 THz (dashed line). 
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